Anomalous density of states in a metallic film in proximity with a superconductor 
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We investigated the local electronic density of states in superconductor-normal metal (Nb-Au) 
bilayers using a very low temperature (60 mK) STM. High resolution tunneling spectra measured 
on the normal metal (Au) surface show a clear proximity effect with an energy gap of reduced 
amplitude compared to the bulk superconductor (Nb) gap. Within this mini-gap, the density of 
states does not reach zero and shows clear sub-gap features. We show that the experimental spectra 
cannot be described with the well-established Usadel equations from the quasi-classical theory. 



At the contact with a superconductor (S), the Andreev 
reflection of the electrons from a Normal metal (N) lo- 
cally modifies the N metal electronic properties, includ- 
ing the local density of states (LDOS) y|. The precise 
LDOS spectra depend on the N-S structure geometry, in 
particular the N metal length, and on the electron trans- 
port regime. In a diffusive N-S junction with a N metal 
shorter than the phase coherence length, one expects a 
fully opened mini-gap, which remains smaller than the 
superconductor's energy gap A 0. In a larger N metal, 
the LDOS shows, within a pseudo-gap, a linear evolu- 
tion with the energy, reaching zero precisely at the Fermi 
level. In the ballistic regime, the same distinction holds 
between a chaotic (mini-gap) and an integrable (pseudo- 
gap) cavity 0. In general terms, a mini-gap shows up 
if every electronic state at the Fermi level can couple to 
the S interface while maintaining quantum phase coher- 
ence. The order of magnitude of the energy gap E g is 
Ti/tar, where tar is the characteristic diffusion time be- 
fore an electron experiences an Andreev reflection. Here, 
we shall consider the case of diffusive N and S metals 
brought in contact through a highly transparent inter- 
face. Then tar. — L? n jD n and the predicted gap E g is 
about the Thouless energy TiD n /L n 2 , where L n and D n 
are the length and the diffusion constant of the normal 
metal. 

Recently, the LDOS of lateral S-N (Nb-Au) structures 
was probed with solid tunnel junctions and (very) low 
temperature STM 0, ||| . These studies focused on the 
pseudo-gap regime in long N metals. Also, a mini-gap 
was observed in a very thin (20 nm) Au layer on top of 
a Nb dot [5j. A good agreement with the quasi-classical 
theory based on the Usadel equations and with the 
Bogoliubov-de Gennes equations Q was obtained. In the 
mini-gap regime, a NbSe2 crystal covered with a varying 
thickness of Au was studied by STM at 2.5 K Q but the 
temperature did not allow the observation of a fully open 
mini-gap. Therefore, the mini-gap regime remains to be 
investigated, in particular its evolution with the normal 
metal size and the possible presence of states within the 
gap. A clear distinction between a fully open gap and a 
pseudo-gap requires a high energy resolution fc^T <C A, 
which can be achieved only at very low temperatures 




FIG. 1: a) Schematic cross section of the full Nb-Au bilayers 
sample, b) STM image (410 x 410 nm 2 ) at 100 mK of the 
sample with a Au thickness of 72 nm. The rms roughness for 
this image is 3.4 nm. 



(T <C IK) and with a large gap A. 

In this paper, we present measurements of the local 
density of states at the N metal surface of S-N bilayers 
with a varying N metal length. We used a very low tem- 
perature STM as the local tunneling probe, which per- 
mitted us to avoid the averaging over the surface inhomo- 
geneities. The achieved high energy resolution enabled us 
to unveil a new phenomenon, namely a non-zero density 
of states appearing in the vicinity of the Fermi level to- 
gether with clear sub-gap structures, and an anomalous 
mini-gap amplitude dependence with the N metal length. 
We demonstrate that these new features cannot be ex- 
plained within the well-established quasi-classical theory. 

We fabricated simultaneously a series of Nb-Au bilay- 
ers, with a fixed Nb thickness and a varying Au thick- 
ness, on a single (6x40 mm 2 ) Si substrate. The Nb 
layer thickness was chosen to be significantly larger than 
the superconducting coherence length in order to avoid 
any effects due to its finite value. For varying the Au 
thickness, another Si wafer was used as a mask and 
moved in situ above the substrate. After depositing the 
120 nm Nb film, the Au film (from 10 to 260 nm) was 
deposited within 15 minutes at a pressure below 10~ 8 
mbar. These conditions minimize the interface contam- 
ination and should preserve the best Nb-Au interface 
transparency. The full Si wafer with the bilayer films 
thus obtained (see Fig. QJi) was cleaved in air to sepa- 
rate the different samples. Individual Au (260 nm) and 
Nb (120 nm) layers were characterized by transport mea- 



TABLE I: Transport properties of the S=Nb (120 nm) and 
N=Au (260 nm) films at 10 K. The RRR is the ratio of the 
room-temperature resistivity with the 10 K residual resistiv- 
ity. 





p (fj£l.cm) 


RRR 


l e (nm) 


D (cm 2 /s) £ (nm) 


Nb 


16.2 


3.7 


5.4 


24.7 23.2 


Au 


2.34 


4.7 


36.1 


169 60.8 



surements which are summarized in Table [I] 

Each bilayer sample was cooled separately in a STM 
working at 60 mK in an upside-down dilution fridge ||. 
A fresh-cut Pt-Ir wire was used as the STM tip. A typical 
STM image (of the 72 nm sample) acquired at very low 
temperature (100 mK) is shown in Fig. We observe a 
polycrystalline structure with a grain size of 30 - 50 nm. 
The same grain size was observed in every sample, except 
in the thinnest one (10 nm). Moreover, this value is con- 
sistent with the measured elastic mean free path l e n in 
the 260 nm sample. Our Au films are thus in the diffusive 
regime, except for the smallest thicknesses where ballistic 
effects may occur. On every sample, we acquired series 
of I(V) tunnel characteristics at a tunnel resistance of 5 
to 12 MQ, and at numerous places in order to check re- 
producibility. The spectra did not depend on the tunnel 
resistance value and were very reproducible over a sample 
surface with minor details (in particular small asymme- 
tries at zero bias) slightly evolving, and also punctual 
defects on a few locations. The I(V) data were numeri- 
cally differentiated to obtain the differential conductance 
dI/dV(V), which gives the LDOS at the energy eV with 
an accuracy limited by the thermal smearing. A selection 
of the tunneling spectra is shown in Fig. The spectra 
are flat at large bias voltages ( | V" | ^ 2.5 meV) and have 
been normalized as 1 there. 

For the bilayer with the smallest Au thickness (10 nm), 
the spectrum qualitatively resembles a BCS spectrum 
with a gap amplitude very close to the expected bulk Nb 
gap : A ~ 1.5 meV. For larger thicknesses, the spectra 
show a mini-gap that reduces in width with increasing 
Au thickness. There seems to be a crossover from one 
type of spectra to another at a Au thickness between 72 
and 100 nm. For small thicknesses, there is a sharp rise 
to a peak in the LDOS at A, and a relatively slow de- 
crease as we go towards zero bias. For large thicknesses, 
there is a small dip at A, the LDOS then rises slowly to 
a peak and falls off rapidly inside the mini-gap. This dip 
at A actually implies that the superconducting gap in 
the bulk of Nb is not significantly influenced by a large 
thickness of Au. 

Inside the mini-gap, the density of states is strikingly 
different from zero and features a few sub-gap structures. 
The observed sub-gap features do not scale with the Nb 
energy gap A, as their energy position evolves with the 
normal metal length L n . At the Fermi level (V = 0), the 
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FIG. 2: Tunneling density of states measured at 60 mK at 
the Au surface of Nb-Au bilayer samples with a varying Au 
thickness L n . Data from the 130 and 200 nm samples are not 
shown for ease of reading. 

mean LDOS monotonously increases with the Au thick- 
ness. This effect is therefore not a ballistic one which 
would be restricted to the smaller thicknesses. 

We checked the energy resolution of our STM by mea- 
suring the LDOS at the surface of plain Nb and Al films 
and fitting the data with a BCS density of states without 
any depairing parameter but with a thermal smearing. 
The fit was very good and gave an effective temperature 
of 210 mK which corresponds to an energy resolution of 
36 /ieV 0. This proven high resolution rules out the 
thermal smearing as a cause of the large (up to 50%) 
zero bias conductance observed here. 

We compared our data to the predictions of the quasi- 
classical theory in the diffusive regime 0, El , by solving 
the Usadel equations with the help of a numerical code 
from W. Belzig et al. Q. The input parameters for this 
calculation are the lengths L n ^ s of N and S metals in units 
of the related coherence lengths £„ iS = ^hD n ^ s /2A 7 the 
specific interface resistance rs, the inelastic 7i„ and spin- 
flip 7s/ scattering rates in the N and S metals, the elec- 
tronic properties mismatch parameter 7 = p s £,s/pn(,n 
(oc ^/l e ,n/le,s)- Here p„, s is the normal-state resistiv- 
ity in N or S. We considered the interface transparency 
as perfect and we took into account a thermal smearing 
with an effective temperature of 210 mK. 

We adjusted the calculation parameters to fit the tun- 
neling spectra. Fig. displays a comparison of exper- 
imental spectra with two sets of calculated curves. We 
first assumed that the inelastic and spin-flip scattering 
lengths are much larger than the Au layer thickness, so 
that we can neglect these processes (Fig. dotted lines). 
In order to recover the observed LDOS peak amplitude 
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FIG. 3: Comparison of the experimental data for L n — 10 ; 50 
; 100 ; 160 and 260 nm (full line) with the calculated spectra. 
The calculation parameters are L s /£ s = 5.17, A = 1.57 meV, 
T = 210 mK, 7 = 0.6, r B = 0, and £ n = 60.8 nm. The dotted 
line curve stands for -yi n = and the dashed line curve for 7i n 
= 0.05 A/h. The curves have been evenly shifted for clarity. 



and position, we had to assume a mismatch parameter 7 
value of 0.6 instead of 2.6 as estimated from the trans- 
port measurements, keeping the other parameters match- 
ing precisely the measured values. This discrepancy may 
be related to the Nb-Au interface roughness which is ex- 
pected to affect the Andreev reflection rate, and thus 
the induced superconductivity in the Au layer The 
overall spectral shapes are qualitatively reproduced, ex- 
cept for the peak of the 50 nm spectra which is notice- 
ably different from the experimental data. At large L„, 
the predicted mini-gap is clearly smaller than the one 
observed experimentally. Moreover, we always obtain a 
fully opened gap with a zero LDOS at the Fermi level. 
This is in clear disagreement with our experimental re- 
sults, where we always get a non-zero conductance at 
zero bias. The dependence of the mean free path with 
the Au layer thickness L n was taken into account 
through the related variation of the characteristic length 
and of the mismatch parameter 7. This modified only 
slightly the calculated spectra in the regime L n < 100 nm 
and did not improve the fit. 

With the hope to obtain a better fit, we included a non- 
zero inelastic scattering rate in the calculation parame- 
ters (Fig. |3]dashed lines). We achieved a good agreement 
of the calculated LDOS with the measured one, but only 
at the precise Fermi level, by choosing j in = 0.05A/?i. 
This value corresponds to an inelastic mean free path of 
about 370 nm, which is significantly smaller than the ex- 
pected value of about 2 fim. At large thicknesses, the 
calculated mini-gap remains smaller than in the experi- 




FIG. 4: Variation of the quantity (E g /A)~ 1/2 with the N 
metal length L n ; the mini-gap E g being defined by the inflec- 
tion point in dl /dV(V). The dotted line shows the mini-gap 
calculated from Fig. [3] parameters (7^ = 0). 



ment. Fitting the mini-gap width requires a significant 
modification of the calculation parameters such that the 
calculated LDOS close to the Fermi level also diminishes. 
This makes the disagreement between the experimental 
data and calculated LDOS even stronger for a fixed value 
of 7i„. We also tried to include a spin-flip scattering 0,0 
or a small interface resistance in the calculation, but it 
did not improve the fit. 

From the experimental tunneling spectra, we extracted 
the mini-gap E g as the energy at the inflection point, 
i.e., with a maximum slope d 2 I/d 2 V. The quantity 
(Eg/A)- 1 ' 2 is plotted in Fig. H as a function of the N 
metal length L n , so that a Thouless energy scaling should 
appear as a linear behavior. The calculated mini-gap is 
close to this scaling : E g = 1.3xl0 4 /(76 + L n ) 2 meV, 
L n being in nm units. As naively expected, the numer- 
ator value is of the order of TiD n — 1.09xl0 4 meV.nm 2 , 
based on Table 1 data. The extra length appearing in 
the numerator expresses the fact that the Andreev re- 
flection happens in the superconductor over this length 
scale. It also makes E g extrapolate to A as L n goes to 
zero. The experimental data follows this behavior in the 
small thickness regime L n < 100 nm. At large L n , the 
discrepancy between theory and experiment is striking. 

Up to now, we described samples with an exceptional 
interface transparency. We have also fabricated similar 
Nb-Au samples but with a degraded interface, thanks 
to an Ar etch of the Nb surface performed just before 
a 20 nm Au deposition. Different Ar etch times were 
used, keeping Ar pressure (3.10 _4 mbar) and beam volt- 
age (100 eV) fixed. Compared to Fig. |3data at similar 
thickness, Fig. ^tunneling spectra show a drastic change 
in the LDOS. A mini-gap is observed, with a decreasing 
amplitude as the etching time is increased. We inter- 
pret this behavior as due to a decrease of the interface 
transparency induced by the surface milling. This en- 
hances the contribution of specular reflection of electrons 
as compared to the Andreev reflection, which makes the 
time tar increase and the mini-gap E g decrease. Inter- 
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FIG. 5: Spectra acquired at the surface of Nb-Au bilayers 
with fixed Nb and Au lengths L a = 60 nm, L n = 20 nm and 
an Ar etched Nb-Au interface. From large to small mini-gap, 
the etch times are 30 s, 1 min. 30 s and 3 min. Inset : fit 
of the 3 min. etched sample data with the Usadel equations 
solution with L s /£ s = 2.5, L„/£« = 0.7, A = 1.47 meV, T = 
210 mK, tb = 6 mQ.cm 2 , 7 = 1.1, ji n — 0. 

estingly, these spectra also show a non-zero density of 
states at the Fermi level, which seems insensitive to the 
variation of the interface transparency. Although we were 
able to qualitatively fit the mini-gap width and the peak 
profile with the Usadel solutions by taking into account 
a non-zero interface resistance (see Fig. [S] inset), it was 
again impossible to fit the whole spectra, especially the 
non-zero LDOS near the Fermi level. A moderate inter- 
face transparency thus causes an anomalous LDOS, in a 
similar way as a large N metal length. 

The observed LDOS, especially its amplitude at the 
Fermi level and the mini-gap width evolution, cannot be 
described by the quasi-classical theory, even after assum- 
ing an extremely strong phase-breaking scattering. In 
the diffusive regime of relevance here, neither the angular 
dependence of the Andreev reflection at an imperfectly 
transparent N-S interface 01 nor the presence of surface 
electron states 14] should play a role. In terms of elec- 
tron trajectories, a finite LDOS at the Fermi level means 
that a significant fraction of the electron states present 
at the Au surface do not Andreev-reflect at the Nb-Au 
interface on the time scale of the phase coherence time. 

In our opinion, this new effect was not observed before 
because previous studies either did not feature the rele- 
vant thickness range and a confined geometry where the 
mini-gap enables a clear observation of additional low- 
energy single particle states 0, @ > or were carried out at 
a higher temperature with insufficient energy resolution 
in the measured LDOS [8j. Still, our observations are 
somewhat reminiscent of some unexplained features re- 
ported on the LDOS of presumably quite disordered 20 
nm N metal ridges lying over a Nb surface jfj. There 
also, a clear filling of the LDOS at the Fermi level was 
observed in the absence of any mini-gap width evolution. 
Let us point out that both this and our systems are far 



from the strongly disordered regime where a non-zero 
LDOS for all energies is predicted 0|. Rather, our films 
have the typical polycrystalline structure of a weakly dis- 
ordered thin film. The elastic diffusion in the N metal is 
not homogenous as in an amorphous metal but controlled 
by the interfaces between ballistic grains. The possible 
confinement of electron states within a grain can have 
serious effects on the LDOS in the proximity supercon- 
ductivity regime of interest here. A significant fraction of 
electron states may be sufficiently decoupled from the N- 
S interface to avoid any Andreev reflection. These states 
would thus remain insensitive to the superconductivity 
in the S layer and contribute fully to the LDOS. 

In conclusion, we investigated with a high energy res- 
olution the LDOS at the N metal surface of diffusive S-N 
(Nb-Au) bilayers and compared to the theoretical predic- 
tions. We uncovered an anomalous mini-gap width evo- 
lution and a non-zero LDOS near the Fermi level that 
increases as the N metal thickness is increased. The lat- 
ter behavior is also observed at small thickness when the 
interface transparency is reduced on purpose. This new 
effect cannot be described within the quasi-classical the- 
ory. We suggest that the granular structure of the Au 
layer enables the decoupling of some electron states from 
the Andreev reflection at the interface, which results in 
a non-zero LDOS inside the mini-gap. 
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